Reflection of electrons from a domain wall in magnetic nanojunctions 
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Electronic transport through thin and laterally constrained domain walls in ferromagnetic nano- 
junctions is analyzed theoretically. The description is formulated in the basis of scattering states. 
The resistance of the domain wall is calculated in the regime of strong electron reflection from the 
wall. It is shown that the corresponding magnetoresistance can be large, which is in a qualitative 
agreement with recent experimental observations. We also calculate the spin current flowing through 
the wall and the spin polarization of electron gas due to reflections from the domain wall. 

PACS numbers: 75.60.Ch,75.70.Cn,75.70.Pa 



C3 



O 

o 



> 
in 
\o 

o 

CO 

o 

-I— > 
c3 



i 

C 

O 

o 



X 



I. INTRODUCTION 

There is a growing interest in the resistance and mag- 
netoresistance associated with domain walls (DWs) in 
metallic ferromagnets. 1 Owing to recent progress in nan- 
otechnology, it became possible now to extract a single 
DW contribution to electrical resistancepi*2iii4 Surpris- 
ingly, it turned out that the resistance of a system with 
DWs in some cases was smaller than in the absence of 
DWsj2i£ whereas in other cases it was larger This in- 
triguing observation led to considerable theoretical inter- 
est in electronic transport through DWsA^iiSiiiii^ The 
interest is additionally stimulated by possible applica- 
tions of the associated magnetoresistance in magnetoelec- 
tronics devices. 

In a series of experiments the magnetoresistance asso- 
ciated with DWs was found to be very large ^Sii 3 . More- 
over, recent experiments on Ni micro junctions showed 
that constrained DWs at the contact between ferromag- 
netic wires produce an unexpectedly large contribution to 
electrical resistivity, and consequently lead to a huge neg- 
ative magnetoresistance^ It was shown theoretically*^ 
that DWs in magnetic micro junctions can be very sharp, 
with the characteristic width L being of atomic scale. 
This is much less than typical DW width in bulk mate- 
rials or thin films. 

Theoretical descriptions of the transport proper- 
ties of DWs are mainly restricted to very smooth 
DWj2ii2ii£ii£*i2iiii which is more appropriate for bulk fer- 
romagnets. Electron scattering from DWs is then rather 
small and the spin of an electron propagating across the 
wall follows the magnetization direction almost adiabat- 
ically. The additional resistance calculated in the semi- 
classical approximation can be either positive or negative 
(depending on material parameters) and rather small. 
The validity condition for the semiclassical approxima- 
tion is kp^n\L S> 1, where fcpf and kp\, are the Fermi 
wavevectors for the majority and minority electrons, re- 
spectively. 



For sharp DWs, however, scattering of electrons from 
the wall is significant and the semiclassical approxima- 
tion is no longer applicable. Some numerical calculations 
of the magnetoresistance in magnetic nanojunctions have 
been presented in Ref. [20] in the context of the con- 
ductance quantization in micro junctions due to lateral 
confinement. The problem of large magnetoresistance in 
magnetic junctions was also analyzed recently by Tagirov 
et aZ;2i where DW was approximated by a potential bar- 
rier independent of the electron spin orientation. The 
ballistic regime of electron transport through the do- 
main wall has been also considered using some numerical 
simulations^ and ab initio calculationsiSi 2 ^ 2 ^ 

In this paper we consider the case of a thin DW, when 
the condition fcp-fmZ < 1 is fulfilled (the semiclassi- 
cal approximation is not applicable). In the limit of 
kp^n\L <C 1, we formulate the problem as a transmission 
of electrons through a potential barrier. Such a formula- 
tion can be treated analytically. In addition, we restrict 
our considerations to the case of DWs with very small 
lateral dimensions, when only a single quantum channel 
takes part in electronic transport. We show, that the 
magnetoresistance associated with DWs can be rather 
large - up to 70%, depending on the polarization of elec- 
trons. 

In section 2 we describe the model and introduce the 
basis of scattering states. Conductance of a domain wall 
is calculated in Section 3. Spin current flowing through 
DW and spin polarization of the electron gas due to re- 
flections from the wall are calculated in Sections 4 and 
5, respectively. Summary and final remarks are in Sec- 
tion 6. 



II. MODEL AND SCATTERING STATES 

Let us consider conduction electrons described by a 
parabolic band, which propagate in a spatially nonuni- 
form magnetization M(r). The system is then described 
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by the following Hamiltonian: 



h 2 d 2 



(1) 



where J is the exchange integral and er = [a x , a v , a z ) 
are the Pauli matrices. For a domain wall with 
its center localized at z — we assume M(z) = 
[Mo sin (p(z), 0, Mo cos <p(z)], where ip(z) varies from zero 
to 7r for z changing from z = — oo to z = +oo. Let the 
characteristic length scale of this change be L (refereed 
to in the following as the DW width). 

When DW is laterally constrained, the number of 
quantum transport channels can be reduced to a small 
number. In the extreme case only a single conduction 
channel is active. In such a case, one can restrict con- 
siderations to the corresponding one-dimensional model, 
and rewrite the Hamiltonian (1) as 



H 



2m dz 2 



JM z (z)a z - JM x (z)a x 



(2) 



Although this model describes only a one-channel quan- 
tum wire, it is sufficient to account qualitatively for some 
of the recent observations. Apart from this, it can be 
rather easily generalized to the case of a wire with a few 
conduction channels. 

In the following description we use the basis of scat- 
tering states. The asymptotic form of such states (taken 
sufficiently far from DW) can be written as 



XRl(z) 



e ik ^ z +r m e-^ z 

J c -rk v z ) . z < ~L, 



ik | z 



(3) 



4-i „ik^z 
L R] e 



where fc T(i) = [2m(E ± M)] 1/2 /%, with M = JM and 
E denoting the electron energy. The scattering state (3) 
describes the electron wave in the spin majority chan- 
nel incident from z = — oo, which is partially reflected 
into the spin-majority and spin-minority channels, and 
also partially transmitted into these two channels. The 
coefficients and t R ^ are the transmission amplitudes 
without and with spin reversal, respectively, whereas tr-\ 
and r R ^ are the relevant reflection amplitudes. It is worth 
to note that transmission from the spin-majority channel 
at z < to the spin-majority channel at z > requires 
spin reversal. The scattering states corresponding to the 
electron wave incident from z = — oo in the spin- minority 
channel have a similar form. Also similar form have the 
scattering states describing electron waves incident from 
the right to left. 

In a general case the transmission and reflection co- 
efficients are calculated numerically, as described in the 
next section. When kp^^L <C 1, then the coefficients 
can be calculated analytically. Upon integrating the 
Schrodinger equation Hip = Etp (with the Hamiltonian 



z > L, 



given by Eq. (2)) from z = - 
L <C 5 <C one obtains 



5 to z = + S, and assuming 




z=+5 



dz 



\a xX (z = 0) = (4) 



z = -8. 



for each of the scattering states (for clarity of notation 
the index of the scattering states is omitted here), where 



dz M x (z). 



(5) 



Equation (4) has the form of a spin-dependent condition 
for electron transmission through a <5-like potential bar- 
rier located at z = 0. To obtain this equation we also 
used the condition k^^L <C 1, which is opposite to the 
condition used in the semiclassical approximation. The 
magnitude of the parameter A in Eq. (5) can be estimated 
as A ~ JM L/h = ML/U. 

Using the full set of scattering states and the condition 
(4) , together with the wave function continuity condition, 
one finds the transmission amplitudes 



t 



RHl) 



= t 



i|(T) - 



t 



f 



AT(l) "I4(T) 



( n +^)2 + 4A 2! 
(«T +"|) 2 +4A 2 



(6) 



where v-m) = ^fc-f(l)/ TO denotes the electron velocity in 
the spin-majority (spin-minority) channel. 

According to Eq. (6), the magnitude of spin-flip trans- 
mission coefficient can be estimated as (for simplicity we 
omit here the state indices) 



It' I 



Xv 



v 2 + \ 2 



Men 



sf so 



M 2 



(k F L) 2 , (7) 



where Ef = h 2 k F /2m and £o = h 2 /mL 2 . For kpL <C 1 
one finds £q £f- Thus, taking ep ~ M, one obtains 



< 1. 



(8) 



Accordingly, a sharp domain wall can be considered as 
an effective barrier for the spin-flip transmission. On the 
other hand, the probability of spin conserving transmis- 
sion is much larger, \t/t^\ ~ EfSq/M 2 ^ 1. This means 
that electron spin does not follow adiabatically the mag- 
netization direction when it propagates through the wall, 
but its orientation is rather fixed. 

It is worth to note, that the conservation of flow in 
the spin-dependent case considered here has the following 
form 



(i-kml 2 ) 



f 

r R-\ 



v l \tR]\ 2 +«T tm ' ( 9 ) 



and also analogous equations for the other scattering 
states. 



3 



III. RESISTANCE OF THE DOMAIN WALL 



To calculate conductance of the system under consid- 
eration, let us start with the current operator 



j(z) = e^(z)vip(z), 



(10) 



where v is the velocity operator, whereas ip^(z) and tp(z) 
are the electron field operators taken in the spinor form. 
Accordingly, the form of Eq. (10) implies summation 
over spin components. Using the expansion of tp( z ) over 
the scattering states (3) and carrying out the quantum- 
mechanical averaging, one obtains the following formula 
for the current 



3{z) 



dk 
2^ 



de 
2^ 



e fe " G n (k,s) xUz)vXn(z) 



(11) 

where n is the index of scattering states (n = 
R |, R i, L |, and L J.) and i] = + . The matrix ele- 
ments of the velocity operator v = —(ih/m) d/dz in the 
basis of scattering states have the form 



vr-\ = (R T | v | R T) = vi \t m \ 2 + w T 



v Rl = (Rl\v\Rl) = n 



r\ 



v L] = {L]\v\L]) = -w T \t L] | - v i 

. 2 



v Ll = (L I | v | L I) 



-V-\ 



Li 



(12) 



Finally, the retarded Green function G n (k,s) in Eq. (11) 
is diagonal in the basis of scattering states. 

When the transmission of electrons through the bar- 
rier is small, one can assume that the chemical potential 
drops at the wall and is constant elsewhere, /i = fj, R for 
z < and fi = /i£ for z > 0. This corresponds to the 
voltage drop U = (fi R — ^i)/e across the domain wall, 
whereas the resistance of the wire parts outside the wall 
can be neglected. The Green function G R ^(k, e) acquires 
then the following simple 1 form 



G R ^(k,e) 



1 



(13) 



where £R-\(k) = h 2 k 2 /2m— M. The other components 
of the Green function have a similar form. 
After integrating over e in Eq. (11) we obtain 

j{z) = e j ^{x R] {z)vxm{z) 0[hr- £R T (fc)] 

+X-RI ( z ) * XRl ( z ) [VR - £Rl ( k )} 
+X L ^ z ) i, XL-\{z) 9[^ L -e iT (fc)] 

+x{ l (z)i X Li(z)6[n L -E Ll (k)}}. (14) 

Since the current does not depend on z due to the charge 
conservation law, it can be calculated at arbitrary point, 
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FIG. 1: Relative conductance of the domain wall as a function 
of its width L in a magnetic quantum wire, calculated for 
indicated values of the parameter p = M/ef- 



say at z = 0. Apart from this, the contribution from the 
states with e R ^i)(k), ££jQ)(fc) < nrin(/i£, /j,r) vanishes 
and only the states in the energy range from min(/i£, fi R ) 
to ma,x(fiL, Hr) contribute to the current. 

Using Eqs. (3) and (12) to (14), one obtains the con- 
ductance G as a linear response to small perturbation 
(t/^0), 



G = 



2irh 



R"i 



Rl 



(15) 

where all the velocities and transmission coefficients are 
taken at the Fermi level. 

When kpK^L <C 1, then taking into account Eq. (6), 
one can write the conductance in the form 



G = 



4e 2 Vfi^ (w T + Vl f + 2A 2 



uf + 



-I) 



( U T + v lY 



4A ; 



(16) 



In the limit of v\ = «| and A — > 0, we obtain the conduc- 
tance of a one-channel spin-degenerate wire, Go = e 2 /irh. 
In the regime of ballistic transport Go is also the conduc- 
tance of the investigated system without DW. 

Variation of the conductance G with the wall width 
L (Fig.l) was calculated from Eq. (15), with the trans- 
mission coefficients determined numerically. Thus, the 
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main wall with finite (nonzero) thickness. This follows 
directly from the weak minimum in some of the curves in 
Fig. 1 (see also Fig. 2). The existence of this minimum is 
related to the sign of the second derivative of the function 
G(A) in Eq. (16), calculated at A = (the first derivative 
vanishes there). In our simple model, the corresponding 
sign is negative for (uj, /«•[■) < {vjv^) cr = 2 — \/3 — 0.268 
(or, equivalently, p < p c ~ 0.866), and positive for 
( u i/ u t) < { v i/ v '[)cr (i-e., p > p c )- When A increases 
from A = 0, the conductance decreases in the former 
case and increases in the latter one. On the other hand, 
we know that for thick domain walls the conductance 
increases with increasing wall thickness. Thus, the mini- 
mum should occur for the curves corresponding to p < p c . 
In the case of strong polarization, p > p Cl the main con- 
tribution to the conductance is associated with the spin- 
flip transmission through the domain wall, and the con- 
ductance increases monotonously with the width of the 
domain wall, in accordance with Eq. (8). 



IV. SPIN CURRENT 



FIG. 2: Comparison of the results obtained from direct nu- 
merical calculations (solid line) and and from Eq. (16) valid 
at k F L < 1 (dotted line) for p = 0.7. 



results shown in Fig. 1, are valid for arbitrary value of 
kpL. The numerical modeling has been done by direct 
calculation of the spinor wave function using Eq. (2), 
starting at z > L in a form of two transmitted spin up 
and down waves with arbitrary numerical coefficients. 
Then we restored the function in the region z <C — L 
and, by numerical projecting the obtained spin compo- 
nents on the right- and left-moving waves (in accordance 
with Eq. (3)), we found the amplitudes of incident and 
reflected waves. 

In the limit kpL <C 1, the results shown in Fig.l 
should coincide with those obtained from the formula 
(16). Comparison of the results obtained from direct nu- 
merical calculations and those obtained from Eq.(16) is 
shown in Fig. 2. Indeed the results coincide for kpL <C 1, 
whereas at larger values oikpL the deviations are large 
and grow with increasing kpL. 

The conductance in the presence of a domain wall is 
substantially smaller than in the absence of the wall. Ac- 
cordingly, the associated magnetoresistance can be large. 
For example, for p=0.9 in Fig.l the magnetoresistance is 
equal to about 70% (which corresponds to G/Go = 0.6). 

It should be noted that in a real magnetoresistance 
experiment on magnetic semiconductor nanowires, for 
which the inequality kpL 3> 1 can be easily fulfilled, one 
can have more than one domain walls. Accordingly, the 
magnetoresistance effect can be significantly enhanced. 

It is also worth to note that the resistance of an abrupt 
domain wall can be smaller than the resistance of a do- 



When the electric current is spin polarized and when 
there is some asymmetry between the two spin channels, 
the flow of charge is accompanied by a flow of spin (an- 
gular momentum). The z-componcnt of the spin current 
can be calculated from the following definition of the cor- 
responding spin-current operator 



J z (z) = tp t (z)a z vif)(z), 



(17) 



which leads to the following average value 

J z (z) = -z£ ^ e ^G„(fc,s) X l(z)a z v Xn (z). 

(18) 

After carrying out the calculations similar to those de- 
scribed in the preceding section, one arrives in the linear 
response regime (limit of small bias voltage U) at the 
following formulas for the spin current J z : 



J z {z < -L) 



J z {z > L) 

Using Eqs. (6) we find 



'R\ 



J z (z>L) = - — 



("T+ W l) + 4A2 



(19) 



• (20) 



(21) 



and J z {z < —L) = —J z (z > L). The magnetic torque 
due to spin transfer to the magnetic system within the 
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FIG. 3: Relative spin conductance of the wire with a domain 
wall calculated as a function of its width L for indicated values 
of the parameter p. 



domain wall is determined by the non-conserved spin cur- 
rent 



T{U) = 



16eU 



A 2 {v* 



(«T +«l) 2 + 4A s 



(22) 



It should be noted that spin-flip scattering due to DW 
does not allow to separate spin channels like it was in the 
case for homogeneous ferromagnets. If we define now the 
spin conductance G s as G s — J z /U, then one can write 
for z > 



G s = 



e 
itTi 



x- 



0? - w l) 



(« T + n ) 2 + 4A2 



(23) 



Thus, G s is negative for z > and positive for z < 0. 

In a nonmagnetic case we have v-\ = and therefore 
G s = 0. In the case considered here, G s = when there 
is no DW. Let us introduce the spin conductance for one 
(spin- up) channel only, G s q — e/2nh. The relative spin 
conductance in the presence of DW, G s /G s o, calculated 
using Eq. (19) and with numerically found transmission 
coefficients, is shown in Fig. 3 as a function of the DW 
width L and for the indicated values of the parameter p. 



It corresponds to the spin current outside the region of 
the domain wall. The spin current inside the wall is not 
conserved because of the spin-flip transitions. 

In accordance with Eqs. (19), (20) and (6), the nonzero 
spin current in a one-channel wire with domain wall is 
due to a difference in spin-flip transmissions for spin-up 
and spin-down channels: the corresponding transmission 
coefficient turns out to be larger for faster (majority) 
electrons. 



V. SPIN POLARIZATION DUE TO DOMAIN 
WALL 

Spin dependent reflections from the wall lead to ad- 
ditional spin polarization of the system near the wall. 
The distribution of spin density created by the wall can 
be calculated using the basis of scattering states. The 
z-component of the spin density in the equilibrium situ- 
ation (U — 0) is 



S z (z) 



dk 
2^ 



de 
2^ 



e 4E "G„(fc,e) xL(z)a, X n(z). 

(24) 



The above formula contains a constant part correspond- 
ing to the spin density in the absence of DW, as well as 
the z-dependent part 8S z {z) created by the wall, 

6S z (z) = — I dk cos(2fc|z) 

dk cos(2fc|z), (z < —L), 

hi 

r-kp 



— / dk r L t cos(2fctz) 

7T JO 



dk Tl[ cos(2fc|z), (z > L). (25) 



The dependence of the spin density S z on the distance 
from the wall is shown in Fig. 4. The spin dependent 
reflections from the wall create spatial oscillations of the 
electron spin density. These oscillations are similar to 
the Friedel oscillations of charge in a nonmagnetic metal. 
However, one should point out here that in addition to 
the above calculated spin polarization, there is also a 
nonequilibrium spin polarization due to flowing current^ 



VI. SUMMARY AND CONCLUDING 
REMARKS 

We have presented in this paper a theoretical descrip- 
tion of the resistance of a magnetic microjunction with 
a constrained domain wall at the contact. In the limit of 
kpm^L <C 1, the electron transport across the wall was 
treated effectively as electron tunneling through a spin- 
dependent potential barrier. For such narrow and con- 
strained domain walls the electron spin does not follow 
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channel. Accordingly, the system was described by a one- 
dimensional model. However, such a simple model turned 
out to describe qualitatively rather well the basic physics 
related to electronic transport through constrained do- 
main walls, although the magnetoresistance obtained is 
still smaller than in some experiments. In realistic situ- 
ations one should use a more general model. When the 
domain wall does not cause transition between different 
channels, then the description presented here can be ap- 
plied directly to the multichannel case by simply adding 
contributions from different channels. 

A domain wall leads to spin dependent scattering of 
conduction electrons. Therefore, it also leads to a net 
spin polarization at the wall, which oscillates with the 
distance from the wall, similarly to Friedel oscillations of 
charge density near a nonmagnetic defect in a nonmag- 
netic metal. We have calculated the equilibrium compo- 
nent of this spin polarization. 

It should be also pointed out that our description ne- 
glects electron-electron interaction. Such an interaction 
is known to be important in one-dimensional systems, 
particularly in the limit of zero bias. The interaction 
may lead to some modifications of the results in a very 
small vicinity of U = 0, but we believe that the main 
features of the magnetoresistance will not be drastically 
changed. 



FIG. 4: Distribution of the spin density of electron gas near 
a domain wall for different values of p. 



Acknowledgments 



adiabatically the magnetization direction, but its orien- 
tation is rather fixed. However, the domain wall produces 
some mixing of the spin channels. 

The calculations carried out in the paper were re- 
stricted to a limiting case of a single quantum transport 



We thank P. Bruno for very useful discussions. 
This work is partly supported by Polish State Com- 
mittee for Scientific Research through the Grant 
No. PBZ/KBN/044/P03/2001 and INTAS Grant 
No. 2000-0476. 



A. D. Kent, J. Yu, U. Riidiger, and S. S. P. Parkin, J. 
Phys. Cond. Matter 13, R461 (2001). 

K. Hong and N. Giordano, J. Phys. Condens. Matter 13, 
L401 (1998). 

U. Riidiger, J. Yu, S. Zhang, A. D. Kent, and 

S. S. P. Parkin, Phys. Rev. Lett. 80, 5639 (1998). 

A. D. Kent, U. Riidiger, J. Yu, L. Thomas, and 

S. S. P. Parkin, J. Appl. Phys. 85, 5243 (1999). 

J. F. Gregg, W. Allen, K. Ounadjela, M. Viret, M. Hehn, 

S. M. Thompson, and J. M. D. Coey, Phys. Rev. Lett. 77, 

1580 (1996). 

N. Garcia, M. Munoz, and Y.-W. Zhao, Phys. Rev. Lett. 
82, 2923 (1999). 

U. Ebels, A. Radulescu, Y. Henry, L. Piraux, and 

K. Ounadjela, Phys. Rev. Lett. 84, 983 (2000). 

G. Tatara and H. Fukuyama, Phys. Rev. Lett. 78, 3773 

(1997). 

P. M. Levy and S. Zhang, Phys. Rev. Lett. 79, 5110 (1997). 
R. P. van Gorkom, A. Brataas, and G. E. W. Bauer, Phys. 



Rev. Lett. 83, 4401 (1999). 

11 P. A. E. Jonkers, S. J. Pickering, H. De Raedt, and 

G. Tatara, Phys. Rev. B 60, 15970 (1999). 

12 G. Tatara, Int. J. Mod. Phys. B 15, 321 (2001). 

13 R. Danneau, P. Warin, J.P. Attane, I. Petej, C. Beigne, 
C. Fermon, O. Klein, A. Marty, F. Ott, Y. Samson, and 
M. Viret, Phys. Rev. Lett. 88, 157201 (2002). 

14 H. D. Chopra and S. Z. Hua, Phys. Rev. B 66, 020403(R) 
(2002). 

15 P. Bruno, Phys. Rev. Lett. 83, 2425 (1999). 

16 G. G. Cabrera and L. M. Falicov, Phys. Status Solidi B 
61, 539 (1974); 62, 217 (1974). 

17 A. Brataas, G. Tatara, and G. E. W. Bauer, Phys. Rev. B 
60, 3406 (1999). 

18 E. Simanek, Phys. Rev. B 63, 224412 (2001). 

19 V. K. Dugaev, J. Barnas, A. Lusakowski, and L. A. Turski, 
Phys. Rev. B 65, 224419 (2002). 

H. Imamura, N. Kobayashi, S. Takahashi, and S. Maekawa, 
Phys. Rev. Lett. 84, 1003 (2000). 



7 



L. R. Tagirov, B. P. Vodopyanov, and K. B. Efetov, Phys. 

Rev. B 65, 214419 (2002); Phys. Rev. B 63, 104428 (2001); 

L. R. Tagirov, B. P. Vodopyanov, and B. M. Garipov, J. 

Magn. Magn. Mater. 258-259, 61 (2003). 

J. B. A. N. van Hoof, K. M. Schep, A. Brataas, 

G. E. W. Bauer, and P. J. Kelly, Phys. Rev. B 59, 138 

(1999). 



J. Kudrnovsky, V. Drchal, C. Blaas, P. Weinberger, 
I. Turek, and P. Bruno, Phys. Rev. B 62, 15084 (2000). 
J. Kudrnovsky, V. Drchal, I. Turek, P. Streda, and 
P. Bruno, Surf. Sci. 482-485, 1107 (2001). 
B. Yu. Yavorsky, I. Mertig, A. Ya. Perlov, A. N. Yaresko, 
and V. N. Antonov, Phys. Rev. B 66, 174422 (2002). 



